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Protecting quantum information from the detrimental effects of decoherence and lack of precise 
quantum control is a central challenge that must be overcome if a large robust quantum computer is 
to be constructed. The traditional approach to achieving this is via active quantum error correction 
using fault-tolerant techniques. An alternative to this approach is to engineer strongly interacting 
many-body quantum systems that enact the quantum error correction via the natural dynamics of 
these systems. Here we present a method for achieving this based on the concept of concatenated 
quantum error correcting codes. We define a class of Hamiltonians whose ground states are con- 
QQ ' catenated quantum codes and whose energy landscape naturally causes quantum error correction. 

, We analyze these Hamiltonians for robustness and suggest methods for implementing these highly 

■ unnatural Hamiltonians. 
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I. INTRODUCTION 



It is tempting to consider all physical systems as being information processing devices that manipulate through their 
dynamics the information labeled by the system's different states. This view of computation, however, is severely naive 
because it leaves out consideration of whether it is practically possible to robustly store and fault-tolerantly manipulate 
^ ' the information encoded into the physical system in the presence of noise and imprecise control ^1, 2, 3]. For classical 
computers, we are blessed with technologies, like the transistor and magnetic media, whose physics guarantees robust 
classical computation. For quantum computers operating according to the quantum circuit model under reasonable 
\ error models, the issue of the existence of robust quantum computation with faulty components was successfully 
^ addressed, in theory at least, by the methods of quantum error correction [J, la, l6| and the subsequent development 
' of the threshold theorem for fault-tolerant quantum computation 0, H, lol. [ToTllll. [T3 |. In spite of this result, the 
\^ , requirements for active quantum error correction are quite severe, and currently no experiment has demonstrated 
■ break even for the most basic active versions of these techniques in a natural setting. Taking a cue from the classical 
world, however, one might wonder whether there exist (or whether we can engineer) quantum physical systems that 
\^ • robustly store information without the use of active quantum error correction but using only the natural physics of 
' the physical system. Kitaev was the first to suggest using the physics of many-body quantum systems, in particular 
to use a topological order, to provide protection for quantum information Following Kitaev's proposal numerous 



"••^ . to use a topological order, to provide protectic 

O ■ related ideas [M [H, [H, [S IH, HO, HH IH, |23|, M l25| have been put forward to protect quantum information in a 
^ , similar manner, some using topological orders and others not topological in nature. A main feature of these schemes 
is that they use quantum error correcting codes along with an energy landscape of a many-body quantum system to 
protect quantum information from decoherence. Roughly one is looking for many-body systems whose energy levels 
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' are quantum error correcting codes and at the same time enforce energetic penalties for quantum errors acting locally 
on the system. None of these methods, however, touch base with the standard method for achieving fault-tolerant 
quantum computation which relies on concatenated quantum codes (an exception being the recent work o Vidal and 
Aguado ^26].) Here we present a new method for physically protecting quantum information based on concatenated 
coding. 

Concatenated codes have a long use in proving the threshold theorem for fault-tolerant quantum computing 
[Ij liol . [ill . [27j . The basic idea is to use a hierarchy of encoding into quantum error correcting codes to protect 
quantum information. If encoding quantum information into a quantum error correcting code can lessen the chance 
of this quantum information being destroyed by quantum noise, then encoding the encoded quantum information 
into a quantum error correcting code should lessen this chance even more. This is the basic idea of concatenated 
quantum error correction. By concatenating r times one decreases the probability of quantum information being 
destroyed doubly exponentially in r. This translates into suppressed of quantum faults at a rate that is exponential 
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in the number of qubits used in the concatenated code. Concatenation is highly reminiscent of renormaHzation group 
methods [28*1 . 

Here we present a class of Hamiltonians, that we call quantum concatenated code Hamiltonians, whose ground 
states are concatenated quantum codes and that have significant protective qualities from decoherence due to the 
energy landscape of the Hamiltonian. The first quality that allows the quantum information to be protected from 
decoherence is the fact that the ground state is separated from the first excited state by a constant gap and the 
smallest quantum error that acts nontrivially on the degenerate ground state contains error that act on qubits 
(where d is the distance of the base code and r is the number of levels of concatenation.) We will show that this 
implies that for local perturbations below a certain threshold strength, the ground state degeneracy is split only 
exponentially in the number of qubits in our system, i.e. as e~"" where n is the number of qubits. We will thus show 
that for local perturbations below a threshold strength there is a phase of the system where quantum information 
encoded into the ground state is protected. Our detailed analytic derivation of this property is the main contribution 
of this paper. This protection is similar to the setting in the toric code Hamilonian [l3|. Further, and beyond the 
toric code in two spatial dimensions, the energy landscape of our code has more than just a constant energy barrier 
between the ground state and the first excited state. Instead quantum errors that affect more qubits lead to states 
of ever increasing energy. Thus there is an energy barrier that errors that destroy the information encoded into the 
ground state must overcome in order to destroy information encoded into the system. This is similar to the toric 
code in four-spatial dimensions 17], and provides an example of a self-correcting quantum system [25|. Thus, whereas 
excitations of the toric code Hamiltonian in two-dimensions during computation can disorder information encoded 
into the ground state of this system, such excitations on a quantum concatenated code Hamiltonian will not destroy 
the information. 

Quantum concatenated code Hamiltonians in their raw form are not reasonable Hamiltonians for a naturally oc- 
curing physical system, since they involve many-qubit interactions on the order of the size of the system. Thus these 
Hamiltonians must somehow be engineered into a physical system. Here we outline one possible approach to achieving 
this engineering. This approach involves creating effective Hamiltonians using strong pulsed interactions. We believe 
this is one of many possible approaches toward designing systems with quantum concatenated code Hamiltonians, 
but leave other methods as an open problem for further research. 

In this work we concentrate on the properties of quantum concatenated code Hamiltonians and one very simplistic 
method for constructing these Hamiltonians using one and two qubit interactions. In doing this we leave out some 
vital considerations for the practicality of these methods. The first, and most obvious, is that we do not consider the 
spatial locality in this paper. The second is that our discussion of the robustness of these models is not yet totally 
complete. In particular methods for rigorously showing the robustness of these models under a generic set of noise 
models have not yet been found. Similarly we do not consider fully fault-tolerant constructions, since the methods 
we use to construct quantum concatenated code Hamiltonians are not themselves robust to introducing errors into 
the system. Finally we note that we address only cursorily methods for performing computation on these systems. In 
many ways, then, our results are akin to early results in the field of quantum error correction, before fault-tolerant 
methods were discovered: we have a developed a method for protecting quantum information but do not have a 
full theory of fault-tolerant quantum concatenated code Hamiltonians. Quantum concatenated code Hamiltonians, 
however, are sufficiently powerful in protecting quantum information that we believe they deserve further careful 
study. 

The layout of the paper is as follows. In Section |TT] we introduce the basic idea of concatenated code Hamiltonians. 
We begin this section by briefly reviewing stabilizer quantum error correcting codes and then introduce stabilizer 
code Hamiltonians and their brethren, quantum concatenated code Hamiltonians. In Section IIIII we analyze the 
robustness of the ground state of these quantum concatenated code Hamiltonians with respect to a local perturbing 
interaction. We show that if such perturbations are below a certain threshold strength, then these interactions 
have an exponentially decreasing impact on the ground state as a function of the system size. In Section IIVI we 
briefly discuss how bang-bang simulation techniques can be used to engineer effective interactions corresponding to 
quantum concatenated code Hamiltonians. In Section |V] we briefly discuss how quantum computation on quantum 
concatenated code Hamiltonians can be performed and discuss the self-correcting properties of these Hamiltonians. 
Finally we conclude with a list of the many open problems for quantum concatenated code Hamiltonians in Scction lVll 

II. QUANTUM CONCATENATED CODE HAMILTONIANS 

In this section we will introduce the notion of a concatenated code Hamiltonian. We begin by reviewing stabilizer 
quantum error correcting codes, a particularly nice set of quantum error correcting codes that we will use exclusively 
in our constructions. We then turn to the deflnition of stabilizer Hamiltonians. These are Hamiltonians whose 
ground state is a stabilizer code subspace and which have a significant energetic gap to the first excited state of the 
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Hamiltonian. These Hamiltonians have a long history in quantum computing going back to the work of Kitaev [l3j . 
Our approach is to define a class of stabilizer Hamiltonians related to concatenated codes. We introduce these 
Hamiltonians in this section, deferring until section HVl how to efficiently construct these Hamiltonians using simulation 
methods. 



We will work with stabilizer quantum error correc ting codes. We briefly describe these codes here. A more detailed 
discussion of these codes can be found in references [13 HO, [Ml . 




Let Vi denote the Pauli group on a qubit, i.e. the group represented by Pauli matrices and a possible phase i'^, 
k e {0,1,2,3}: Vi = {i^I,i^X,i^Y,i^Z}. The Paufi group on n qubits, Vn, is just the group represented by tensor 
products of Pauli matrices and a phase , k € {0, 1, 2, 3}. If P = i''Pi ® • • • (g) P„ is a Pauli operator, then the weight 
of this operator is the number of Pi that are not /. On a similar note we call an operator fc-local if it can be written 
as a sum of Pauli operators with weight k. All elements of the Pauli group square to ±/ and all elements of the Pauli 
group either commute with each other, PQ = QP or anti-commute with each other PQ = —QP. 

Stabilizer codes [s^ are defined by the use of a particular subgroup of the Pauli group known as a stabilizer 
group. A stabilizer group is a subgroup of the Pauli group that is abelian and that does not contain — /. Since all 
of the elements of S commute we can simultaneously diagonalize these operators. Let Si, ... , Sn~k denote a minimal 
generating set for a stabilizer subgroup S. Since the stabilizer subgroup does not contain — /, each stabilizer element 
squares to /, and thus has eigenvalues ±1. Simultaneously diagonalizing the Si allows us to label subspaces of the 
Hilbert space of n qubits by the eigenvalues of these operators. In particular, we can define a basis |si, . . . , Sn-k, m) 
where Si € {±1} and m e {0, . . . , 2'' — 1}, such that Si\si, . . . , s„_fe, m) = Si\si, . . . , s„_fe, to). In other words for 
each of the 2""*^ possible eigenvalues of the stabilizer generators we obtain a subspace of dimension 2^^, that we call 
a stabilizer subspace. A stabilizer code is usually defined as the +1 simultaneous eigenspace of the Si operators, 
Sili^c) = IV'c) for \ipc) in this code. Note however, that we could equally as well have chosen any ±1 subspace and 
obtained an equivalent code. 

Given a stabilizer S, the normalizer of S in the Pauli group, J\f, is the set of elements of the Pauli group that 
commute with the stabilizer S. Note that the normalizer M contains the stabilizer S. If we examine the group Af/S, 
this group acts nontrivially on the different stabilizer subspaces. In particular it is isomorphic to a Pauli group on k 
qubits. In other words it is possible to choose these operators in such a way that they act as logical Pauli operators 
on the k qubits encoded into the stabilizer code. We will denote the ith encoded Pauli operator on the logical qubits 
by Pi where P is the corresponding Pauli operator. Note that there is no canonical chose for this labeling, just 
as there is no canonical choice for the generators of the stabilizer group. We can then label the degeneracy of the 
|si, . . . , s„_fe, to) basis via logical qubits, |si, . . . , s„_fc, 6i, . . . , bk), where the corresponding encoded Pauli operator 
acts accordingly, Pi\si,. . . ,s„_fc,6i, . . . ,6^) = |si, . . . , Sn~k) ^E) |6i) (8) • • • Pz\bi) (g) • • • (g) \bk). 

We are now in a position to describe the error correcting properties of a stabilizer code. Since the dimension of 
a stabilizer with n — k minimal generators is 2*"', a stabilizer code encodes k qubits into n qubits. Suppose that 
E anticommutes with a stabihzer generator Si. Then Si{E\si, . . . , s„, 6i, . . . , 6„)) — —ESi\si, . . . , s 

n 1 ^1 ; ■ • ■ ; ^n) — 

—SiE\si, . . . , Sn, 6i, . . . , bn), SO we see that the effect of E is to flip the sign of eigenvalue of the stabilizer generators 
with which E anticommutes. This implies that all such E^s are detectable errors: if we start in the all +1 simultaneous 
eigenspace of the stabilizer generator and if E is applied to this state, then we can measure the eigenvalue of the 
stabilizer generators and detect that an error has occurred. Let {Ea} denote a set of Pauli operators that will represent 
errors. Then consider the set made up of pairs of these errors as {EIEi,}. If every E^Ei, is not in the normalizer 
minus the stabilizer, i.e. is a detectable error or in the stabilizer, then we can use the stabilizer code to correct the 
errors {Ea}. We are mostly interested in codes that can correct errors {Ea} made up of all Pauli operators up to a 
weight t. A code that can correct t errors must have a distance 2t + 1 where the distance of the code is the smallest 
weight Pauli operator that acts nontrivially on the code subspace. We label a quantum error correcting code on n 
qubits that encodes k qubits and has a distance d via the standard notation [[n, 

Subsystem codes are a generalization of subs pac e codes in which instead of encoding information into a subspace, 
one encodes this information into a subsystem '25[, [3l|, [s^ [s^, [1^ • For our purposes we will only be concerned with 
subsystem stabilizer codes. A simple way to think about subsystem stabilizer codes is as follows. Take a [[n, k, d]] 
stabilizer code. Now assign r of the logical qubits as gauge qubits. We can then encode quantum information into the 
remaining k — r encoded qubits. These fc — r qubits are a subsystem of the k encoded qubits. If we stipulate that we do 
not care what happens to quantum information encoded into the gauge qubits, then we are encoding into a subsystem 
code, where we care not about preserving information encoded into a subspace, but information encoded into a 
subsystem. Such codes can always be turned into subspace codes, but their recovery operations can be significantly 
simpler than the corresponding subspace codes [2^, [Sg] , since one does not care about errors accumulating on the 
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gauge qubits. We will denote a stabilizer subsystem code with r guage qubits, encoding k qubits, and have a distance 
d by the notation [[n, k, r, d]]. 

We now turn to describing briefly how stabilizer codes are used to obtain thresholds in quantum error correction. 
It is this scheme that we seek to mimic via physical means. Suppose we have a quantum error correcting stabilizer 
code C that encodes a single qubit into n qubits and has a distance d and hence can correct t = [^y^J quantum 
errors. If single qubit errors occur probabilistically and independently with probability p on this code, then quantum 
error correction on this code will fail when greater than t errors occur (assume perfect encoding and decoding for 
now as we are simply trying to give the jist of thresholds in error correction.) This implies that the probability of 
failure is pfaii < Cp*+^ for some fixed constant C. Concatenating a quantum error correcting code with itself is the 
process of using encoded codewords to form a new code. Concatenating a code with itself produces a new code with 
qubits, a distance d^, and a probability of failure less than C(Cp*"'"^)*"'"^. The procedure for correcting an error on 
a concatenated code is to correct the errors on the first level of encoding and then to correct the errors that remain 
on the next level of concatenation. Concatenating r times, uses qubits and has a probability of failure given by 
recursing the above failure probability formula, 

— = ( — ) > where = — . (1) 

P* \P* J C 

Thus if p < (below a threshold) the probability of failing decreases doubly exponentially in the number of levels of 
concatenation and exponentially in the number of qubits. This scaling is the basic principle behind the threshold for 
fault-tolerant quantum computing. We seek to mimic this concatenation method, but now using the energy structure 
of a Hamiltonian. 



B. Stabilizer Code Hamiltonians 



We now turn the construction of Hamiltonians for many qubit systems that have energy levels that are error 
correcting code subspaces and whose energy landscape enforces energetic costs for quantum errors. We begin our 
discussion with what we term stabilizer code Hamiltonians. 

From a stabilizer [[n, code we can easily construct a simple Hamiltonian whose ground state is degenerate 
with this degeneracy corresponding to the stabilizer code space. One way to achieve this (dating back to the work 
described in see also [2^) is to make the Hamiltonian a sum of the stabilizer generators Si, 

Hstah = ""^ X] 

i 

or, if we express this in the |si, . . . , Sn-k, bi, . . . , bk), basis 

-ffstab = ^ ^ ^ Sj|si, . . . , S„_fc, 5i, . . . , 6fc)(i, . . . , S„_fe, 6i, . . . , 6fc|. (3) 

si,....s„_fce{±l}6i,...,6fce{0,l} i 

The eigenvalues of this Hamiltonian are between — -^(n — 1) and — -^(n — 1) in increasing steps of size J. Every energy 
level is at least 2''-fold degenerate and is separated form the nearest energy levels by J. In fact if we label the energies 
of the Hamiltonian by an index i from to n — 1, such that the energy of the ith energy level is Ei — —^{n—l) + Ji, 
then the degeneracy of the zth energy level is gi — 2("T^). Note that since the generators of a stabilizer code are not 
unique, Eq. ^ is not unique. We will always assume a fixed set of generators and logical operators have been chosen 
for a stabilizer code Hamiltonian. Further note that this is definitely not the only manner to take stabilizer elements 
and form a Hamiltonian whose ground state is a stabilizer code. 

Stabilizer Hamiltonians have some remarkable properties that we will now describe. Since we are working with a 
code of distance d, all Pauli operators with weight less than d act as either identity or on the ground subspace. 
Let {E} denote a set of detectable errors, that is elements of the Pauli group which anti-commute with at least one 
stabilizer generator. Then every such error will take elements of the grounds state of the Hamiltonian to an excited 
state, since these operators will always flip the sign of at least one Si when acting on a |si, . . . , Sn-k, bi, . . . , bk) state. 
Thus, for example if our code is a one qubit error correcting code, then the smallest weight Pauli operator that acts 
as a non-trivial operator on the ground state of a stabilzier code Hamiltonian is of weight three. In other words every 
weight two Pauli operator acting on the ground state of this Hamiltonian causes the ground state to be excited to a 
higher energy level. Since our system has an energetic gap of size J every single qubit and two qubit operation "costs" 
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an energy of at least J. In general if we have a [[n, k, d]] code, any weight w < d Pauli operator costs an energy at 
least J. 

The first stabilizer Hamiltonian where the discussion was explicitly cast in terms of quantum error correction was 
Kitaev's toric code Hamiltonian [Ts']. Kitaev's toric code is a stabilizer code where the generators are four qubit 
operators acting on qubits on links of a square lattice with toric boundary conditions. If this lattice is of size , then 
Kitaev's code is a [[2L^,4, L]] code. Recently Jordan et al [23] introduced stabilizer Hamiltonians in order to provide 
energetic protection within the context of adiabatic quantum computing. In particular they considered small [[4, 1, 2]] 
error detecting codes (see also [s^l) as well as the [[5, 1, 3]] single qubit error correcting code stabilizer Hamiltonian. In 
both Kitaev's construction, Jordan et al's constructions, and in the similar constructions of coherence-preserve (also 
known as supercoherent) qubits [IBjIlSl) some forms of decoherence on information encoded into the ground state can 
be suppressed. This is because in all of these models, single, double, or multi-qubit errors excite the ground state to 
an excited state that is of a greater energy by some constant energy gap A. Then if, for example, the bath is a simple 
set of harmonic oscillators at temperature T <C A, these decoherence mechanisms will be exponentially suppressed 
due to the lack of occupancy of bath modes with enough temperature to excite the ground state. 

As an example consider of a stabilizer code Hamiltonian, consider the [|5, 1, 3]] quantum error correcting code. 
This is the smallest stabilizer code that can correct a single qubit error [H, [33|. This code is a nondegenerate code, 
meaning that every error of weight less than 3 is a detectable error and is not an element of the stabilizer. For the 
[[5, 1, 3]] stabilizer code, the stabilizer Hamiltonian can be chosen to be 

H = -'^ {XZZXI + IXZZX + XIXZZ + ZXIXZ) . (4) 

The ground state of this Hamiltonian is the subspace of all +1 eigenvalue eigenstates of the individual terms in this 
Hamiltonian. The logical operators on this stabilizer code may be taken to be P = P®^ (although smaller weight 
choices are possible.) 

Another example of a stabilizer code Hamiltonian, but this time constructed from a subsystem code and not 
a subspace code, can be derived by using the code described in [25]. This is a [[9,1,4,3]] subsystem code which 
encodes a single qubit, uses 4 gauge qubits, and has distance 3. The stabilizer for this code can be generated by 
{XXXXXXIII, IIIXXXXXXX, ZZIZZIZZI, IZZIZZIZZ} and thus the stabilizer Hamiltonian can be chosen 
to be 

= - {XXXXXXIII + IIIXXXXXXX + ZZIZZIZZI + IZZIZZIZZ) . (5) 

The stabilizer subspaces for this code is 2^ — 32 dimensional, representing 5 encoded qubits. However, 4 of these qubits 
are gauge qubits, degrees of freedom which we do not wish to protect from the detrimental effects of decoherence. 
A logical X operation on this subsystem code can be chosen to be XXXIIIIII and a logical Z operation on this 
subsystem code can be chosen to be ZIIZIIZII. Note that in subsystem codes there is degeneracy beyond the logical 
qubit degeneracy and which corresponds to the gauge qubits. 

C. Concatenated Stabilizer Code Hamiltonians 

We now move on to describing concatenated stabilizer code Hamiltonians. The basic idea here will be to construct 
Hamiltonians whose ground states are a concatenated stabilizer code. 

We begin by describing the first concatenation. Assume we are working with a [[n, l,c?]] stabilizer code. Take 
qubits. On blocks of size n of these qubits we can implement a stabilizer code as in Eq. Further, using the logical 
Pauli operators on each of these n stabilizer codes, we can implement a stabilizer code on the encoded information. 

In particular we can define the rth concatenated Hamiltonian recursively 

n 

Hr^J2"r\+HT^b°''°''. (6) 
i=l 

where is the level r — 1 concatenated stabilizer code acting on the ith block of qubits and HXab'^^'^ ^ stabilizer 

code acting on the encoded quantum information of each of these blocks. This Hamiltonian will contain — 1 
terms with Pauli operators as large as weight n. This later fact is the challenge in implementing this Hamiltonian 
physically [40ll4ll. |4^. For now, however, assume that you could implement this Hamiltonian. We will return to the 
question of implementing this Hamiltonian after we analyze properties of this Hamiltonian with respect to quantum 
errors. 
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For concreteness, here is the concatenated Hamihonian at r = 2 for the [[5, 1, 3]] code, where we use the stabihzer 
HamUtonian described in Eq. ^ and the logical operators on this stabilizer code of Pj — P^^; 



2 ^ 

— -jH = X5i+iZ5i+2-^5i+3-'^5i+4 + -'^5j+2-^5i+3^5i+4-'^5j+5 + -'^5i+l-'^5i+3^5i+4^5i+5 + ^5i+l-'^5j+2-'^5j+4^5j+5 

4=0 

+XiX2X3X4X5Z(jZjZsZgZioZiiZi2Zi3Zi4Zi5XiQXijXisXigX20 

+XQXjXsXgXioZiiZi2Zi3Zi4ZirjZieZijZisZigZ2oX2lX22X23X2AX25 

+XiX2X3X4X5XiiXi2Xi3Xi4Xir,ZieZijZisZigZ2oZ2lZ22Z23Z24Z25 

+ ZiZ2Z3Z4Z5X(iX7XsXgXioXiQXnXisXiQX2aZ2oZ2lZ22Z23Z24Z25- (7) 

We have used the notation that Ri denotes the operator R acting on the ith qubit. 

III. ROBUSTNESS OF QUANTUM CONCATENATED CODE HAMILTONIANS WITH RESPECT TO 

LOCAL PERTURBATIONS 

We will now analyze the properties of concatenated code Hamiltonians with respect to serving as a robust quantum 
memory. We will begin by considering the effect of errors on stabilizer Hamiltonians and then use this understanding 
to analyze the concatenated code Hamiltonians. We first discuss these effects on subspace stabilizer code Hamiltonians 
and then discuss how these generalize to subsystem stabilizer code Hamiltonians. 

A. Single Qubit Errors on Distance Three Subspace Stabilizer Hamiltonians 

Consider the effect of single qubit perturbations on a stabilizer subspace code Hamiltonian (not concatenated) for 
a [[n, 1, 3]] code. To begin with consider a perturbation of the form 

TL 

i=l Qe{X,Y.Z} 

where Ri denotes the operator R acting non-trivially only on the ith qubit and x is our perturbation expansion 
parameter. It will be convenient to let maxg^i Xq^i = 1 such that x represents the strength of the strongest single 
qubit term. Note that this implies that \\V\\ < 3n. We will use the degenerate perturbation theory described in 
Appendix [X] This perturbation represents a static single qubit perturbation. We will delay until later discussion of 
more general setups where our qubits are coupled to baths. 

We consider F as a perturbation on i/gtab of Eq. ^ and focus on the effect of this perturbation on both the 
ground state of H and the excited states of iJstab- Recall that these states are made up of stabilizer subspaces, often 
with multiple subspaces corresponding to a highly degenerate eigenvalue. Further we will work under the assumption 
that a; 1 1 1^1 1 ^ ^ so that the eigenvalues are not shifted by an amount comparable to the gap in our unperturbed 
Hamiltonian. 

We will begin by deriving an effective Hamiltonian for the ground state. The ground state of a [[n, 1, 3]] stabilizer 
Hamiltonian is two-fold degenerate with this degeneracy corresponding to the stabilizer code. Let Ho denote the 
projector onto this degenerate ground state. Recall that the smallest weight operator that does not have a vanishing 
action on a stabilizer subspace for a non-degenerate stabilizer code has weight d. 

Begin by considering the first order perturbation, g'^^VG^q'^ — HoVTlg. This term vanishes because V is made 
up of a sum of single qubit operators and each single qubit operator acts to take a stabilizer subspace to another 
stabilizer subspace, since every single qubit error must anticommute with an element of the stabilizer 

x4'^ =xG^°W^^ ^0 (9) 

(We can ignore the case where single qubit errors are degenerate as this corresponds to qubits that are in fix states 
unentangled with the rest of the qubits and can be discarded.) 

Next consider the second order perturbation g[,°Vg'^^ Vg[,"' + G^,^ VgJ," Vg[,"^ + g[," Vg[,"Vg'[,^^ Only the first 
of these can be nonzero because, as we have just said, G^'^VG^'q^ = 0. Explicitly expanding this term wc find that it 
is 

n TT ^ 

noE E Ag,Q.E rig 'e) Y ^ ^B.,R,no. (10) 
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Since the distance d of the the stabilizer codes is greater than 2, every weight 2 Pauh operator has a trivial action 
on stabilizer subspaces, acting either as the identity on the code space or as zero on the code space. Let (T2 denote 
the number of terms in the cross expansion among operators in the above equation that act as identity. Note that 
0-2 is at least 3n, since every Pauli operator squares to identity, but at most 9n^, since there are only 9n^ terms. In 
particular this implies that the second order term is proportional to the projector onto the ground state subspace 
with a constant depending on the specifics of the code being used, 

x'^A^^^ = xHoIio, (11) 

where we can bound bo 

M < y, (12) 

using the fact that the stabilizer Hamiltonian has a energies separated by at least J. Thus we see that the lowest order 
term in our perturbation has the effect of shifting the energy levels by an amount bounded by x'^bo- Since x is the 
strength of the individual qubit couplings, when x -C ^ the energy shift has much less strength than the individual 
qubit couplings. 

If we truncate the above perturbation theory expansion at second order we obtain the effective Hamiltonian for the 
eigenvalues arising from the ground state of 

Hf = EMx) + xXTIo + Mo, (13) 

where Mq is an operator that acts arbitrarily, but whose operator norm is bounded by 

\\Mo\\<j[^)\ (14) 

Note that this implies that this part of the effective Hamiltonian gives rise to an error on information encoded into 
the ground state that is third order in x. It will turn out that there are second order errors that are more important 
than these errors, but already one can see that the effect of encoding is to reduce an error of strength a; to a level 
shift of strength x"^ and an error of strength x^. 

Next consider the effect of xV on the excited states of the stabilizer Hamiltonian. Recall that a given energy level 
for a stabilizer Hamiltonian consists of multiple stabilizer subspaces. Let i denote the ith energy level of the stabilizer 
Hamiltonian and let a label the different stabilizer subspaces that make up this ith energy level. The first order 
perturbation term is given by G\^^VG\^\ This is no longer zero, since a single qubit error can move one between the 
different stabilizer subspaces. Define Qi.a,i3 as an operator with operator norm less than 1 that acts between the ath 
and /3th stabilizer subspaces of the ith energy level. Then 

= xGf Vcf ^ = X ^ ai{a,l3)Qi,^,f3. (15) 

We can bound ai(a,/3) by counting the number of ways that the single qubit errors can take the ath stabilizer 
subspace to the /3th stabilizer subspace. Note that this error acts entirely within the excited state energy manifold 
and does not connect to any other energy levels. It is this property that will render this term largely irrelevant to our 
analysis. 

The second order perturbation term is given by Vcf ^ + Vcf Vcf ^ + Vcf Vcf ^ The first 

of these terms will give rise to a term that acts entirely within the ith energy manifold. Like the first order term it 
can expressed entirely in terms of operators that act between the stabilizer subspaces, 

Vcf Vcf ^ = ^ bi{a, m,a,0, (16) 

a,l3 

where Ri^a,0 is an operator whose norm is bounded by 1 and we can bound bi{a, jS) by the number of ways that 
two single qubit errors (perhaps operating on the same qubit) take the ath stabilizer subspace to the /3th stabilizer 
subspace divided by the energy gap that is at least J (coming from the g|^'' term.) 

Next consider the second and third terms in the second order perturbation term. This is the first term that can 
take an excited energy level and move it to another energy level. Indeed we see that this term can be nonzero when a 
single qubit error takes the ath stabilizer subspace to a /3th stabilizer subspace of the same energy, and then another 
single qubit error takes that /3th stabilizer subspace to another stabilizer subspace of energy different that i. The net 
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effect of this is an operator that takes the ath stabihzer subspace of the zth energy level to the a'th stabilizer subspace 
of the jth energy level. Let Ti j ^^a' denote an operator of maximum norm 1 that acts between the ath stabilizer 
subspace of the «th energy level and the a'th stabilizer subspace of the j ^ ith energy level. Then the perturbative 
terms can be written as 

Gf VcfVci"' + ycf VGp) = E E ^^^■("' (17) 

i^j a.a' 

Note that we can bound bij{a,a') by the number of three qubit processes that cause errors of the form described 
above divided by the energy gap between the ith and jth energy level, which is at least J. Note that this term can 
take the ground state to an excited state of the stabilizer Hamiltonian. Importantly, however, this term is of strength 
on the order of . 

Putting this together we have derived an effective Hamiltonian for the i energy level given by 

Hf = E,n,{x) - x E a«(a, f3)Q^,o.,p + a:' E /5)-^^",/3 + E E + (1^) 

0.^(3 a, (3 i^j a,a' 



where again the error term is bounded by 



l|M.||< J(^)'. (19) 



Recall that the effective Hamiltonian Hf^ is exactly the Hamiltonian whose eigenvectors are the eigenvalues arising 
from the slitting of the zth energy level. Thus summing over all of the effective Hamiltonians will give us exactly 
H = Hq + xV . In particular we have derived 



H = HU,^ + x'bono+J2 



,ct.a.' 

a,/3 i^j a, a' 



+ M, (20) 



where 



\m\<j{^^)\ (21) 

using the fact that each effective Hamiltonian gives rise to perpendicular subspaces and defining 

n-l 

^fstab = E™(^)- (22) 

Note that -ffg^ab ^'^'^^ each new eigenspace in a trivial manner, i.e. only shifts the energy of this energy eigenspace. 
This is not to say that this operator preserves information encoded into the ground state of the original -ffstab, but 
rather that the ground state is changed to a different ground state. In other words, H'^^^^ is not responsible for splitting 
the energy levels. We analyze this term later and show that the new ground state is close enough to the original 
ground state, where close enough will be defined in terms of a quantum error correcting code order parameter [25| . 

Let us focus on how the H wc have derived acts on the ground state. The ground state is acted upon by the first 
two terms and last two terms in Eq. ([^0)1 . For now we will ignore the effects of the first term and focus on the other 
three terms. If our stabilizer code has the basis |si, . . . , s„~i, 6), then we can express a generic operator on our space 
as 

E Wl®L, (23) 

where Wl acts only on the stabilizer tensor component of our states, . . . , Sn-i,b) — {Wl\si, . . . , s„_i)) |6). 

In particular we can do this for the second term and the last two terms in Eq. (j20[) . 



E Vl^L^ x^boTlo + E E ^ 

ij{a,a')Tij^a,a' + M. (24) 

Le{I,X,Y,Z} i^j a, a' 
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When we do this we can further obtain a bound on the strength of Vl' 



f 1 2?7t\ ' 

II^lII <a;2|6o| + :E' maxmax|6o.(0, + ■ (25) 

i^Q a \ J J 

Note that |6oi| is at most ^^j^ by a naive counting of all possible processes that can act as two qubit errors of the 
correct form. Better bounds can be obtained for specific codes using the combinatorics of how errors occur on such 
codes. We perform such an analysis for the five qubit code in Section [III A 21 

It is useful to take a step back and examine what we have shown so far. We have shown that it is possible to describe 
the effect of a single qubit perturbation on a stabilizer Hamiltonian in a form that turns the stabilizer Hamiltonian 
into a new perturbed stabilizer Hamiltonian H^,^^^ and acts on the ground state only nontrivially in order x'^, i.e. of 
strength less than Cx^ for some constant C and some small enough x. Thus if x < ^ this interaction strength is less 
than the single qubit coupling strengths and we have effectively reduced the strength of interactions that can cause 
information encoded into the ground state to be destroyed. Thus we see that by encoding information into the ground 
state of this new Hamiltonian, single qubit perturbations can destroy this information, but do so to a lesser degree 
than if we had not encoded the information. But what about this new perturbed Hamiltonian, H!,^^^7 

We will now show the eigenspaces of -ffgtab actually very close to those of ifstab, in a well defined manner. In 
particular we will calculate the projector onto the ith perturbed subspace, Tli{x) via the perturbation theory describe 
in Appendix [XI To first order, we find that 

U^ix) = n, - a;(Gf Vcf ^ + Vcf ^) + iV„ (26) 

where we can bound 

\m\<2i^-^j . (27) 

Note that the first order term does not vanish. This first (and higher) order terms appearing here imply that the 
eigenstates of -ffg^ab ^'^ longer made up of stabilizer subspaces. What do these eigenstates look like? Suppose we 
take a stabilizer subspace state {ipi) and apply the above projector. Then we will obtain 

n,(x)|V,) = \^,) - xiG'i'WGf'^ + G^^Vcf ^)|V^,) + x^a\^,), (28) 

where we can bound a < 2 (i^) and is some arbitrary state dependent on Ni. Now notice that {G[^'^VGf'^ + 

Gf'''VG[^^)\tpi) consists of single qubit errors acting on the the stabilizer subspace state \ipi). Since our code can 
correct a single qubit error, this component of the wave function represents a correctable error. Further we can bound 
the amplitude of this term by noting that every single qubit error can contribute at most j amplitude to this state. 
In other words, the new energy levels of H^^.^^ are made up of a superposition over the original stabilizer code state 
plus single qubit error terms acting on these states of relative amplitude 0(a;), plus a term whose relative amplitude 
is 0(x2). 

At this point it is important to introduce a quantum error correcting code order parameter [2^, Suppose 
that you have information encoded into an error correcting code and then some noise process occurs on this encoded 
information. For noise processes of a correctable form, i.e. for say single qubit errors, the information encoded into 
our code has not really been destroyed, since there is a measurement we can perform that will allow us to restore the 
quantum information. Thus, given a quantum error correcting code, we can define a quantum error correcting order 
parameter for stabilizer codes as follows. Suppose we are trying to define a quantum error correcting order parameter 
for information encoded into a stabilizer subspace jsi, . . . , s„_i, 5) for a fixed set of Si. Correctable errors on this 
code either act trivially on these states or as detectable errors with particular syndrome measurements. In particular 
for each Pauli error that is correctable and not trivial, E, we can define a syndrome si{E), . . . , s„_i(£') which, if we 
measured these outcomes for the stabilizer generators, would indicate that error E has occurred. Then the quantum 
error correcting order parameter is given by 

D = P(si, . . . , s„_i) + EP{si{E), E{sn-i)), (29) 



where the sum is over all correctable errors. The importance of this error correcting order parameter is as follows. 
Suppose that we encode information into |0) and then an error process occurs that turns this state into a new state 
p. Then {(j)\D^ pD\(/)) has contributions from p that are 1 if the information is correctable, and is otherwise. 
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Consider what this all means for the use of a stabilizer code Hamiltonian as a quantum memory. In the presence of 
a constant and small perturbing interaction of the form we have considered, two effects occur. One is that the order 
X couplings are changed into order couplings. The second effect is that the ground state of the Hamiltonian is 
slightly perturbed. This perturbation produces effects to first order that are correctable, but that to second order are 
not correctable. Thus both detrimental effects have changed first order perturbations to second order perturbations. 



1. Coupling to a Bath 

Finally it is important to note here that while our derivation was for a perturbation from a constant single qubit 
perturbation, our results are actually more generic than this. In particular consider the case where our qubits are 
coupled to a bath. Then our perturbation can be expressed as 

n 

xV = xY^ AQ,,g, (K)Bq^, + i7bath, (30) 

j=l Qi£{X,Y,Z} 

where Bq i and i?bath act on the bath's Hilbert space. Then if we absorb in x the strength of the bath operators into 
Bq^i, such that ||BQ,j|| < 1 and work in the rotating frame of the bath Hamiltonian, then our above derivation will 
carry through. 

Note, however, an important caveat about this procedure. By assuming that we have absorbed the strength of Bq i 
into X we have limited the models for coupling of a quantum system to its environment. In particular if the bath 
operators are not bounded, as for example happens when the bath is a harmonic oscillator, then our procedure will 
not work. This is a restriction in our analysis. We expect, however that it practice it will be less of a problem, since 
in practice our baths do not have initial states that have unbounded energy and that will cause the unboundcdness 
to become a problem. However, a detailed analysis of this problem is beyond our initial investigation of quantum 
concatenated code Hamiltonians. 



2. Example for the Five Qubit Code 

Here we derive the combinatorial factors for the five qubit code [H, [3§| that allow us to explicitly bounded the 
effective logical qubit strength. Recall that Eq. ^ gives a stabilizer Hamiltonian for the [[5,1,3]] subspace code, bo 
is the number of ways that a two qubit error process can produce identity on the stabilizer code. Since the five qubit 
code is a perfect code, the only way that this can happen is when the same two qubit error occurs on a qubit. This 
implies that |6o| < ^■ 

Similarly we can bound boi{a). Recall that the energy levels of a stabilizer Hamiltonian arc between ~^{n — 1) 
and — 1) in steps of n, or Ei — —^{n — 1) + Ji. The number of stabilizer subspaces that make up the ith energy 
level is given by ("7^)- -^Y explicit computation we can see that the largest error process for the relevant two qubit 
errors are ones that cause excitations to the first excited state. In particular we can bound 



maxmax6oi(0, a') 

i^O a' 



< (31) 



since there are 4 stabilizer subspaces to excite to and 3 remaining for the second transition in this second order 
perturbation effect. 

Putting this together we find that the effective strength of our single qubit interactions on the encoded information 
has become 

15 2 2 12 r f^Ox^ ^ 



\\VL\\<-x'+x'- + Ji^—j . (32) 

If we are interested in where this becomes less than x, we obtain a polynomial inequality in terms of 7 = -j 

277 + 60^72<l, (33) 



which is true for 7 < ^fZV^ « 0.0020901. 
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B. Single Qubit Errors on Distance Three Subsystem Stabilizer Hamiltonians 

Above we have consider the effect of single qubit perturbations on a subspace stabihzer Hamiltonian. But what if 
we use a subsystem stabilizer Hamiltonian instead? In this case, errors which act on the gauge qubits do not have 
any effect on information encoded into the logical subsystem. Much of the analysis which we have performed above 
carries over to subsystem stabilizer code Hamiltonians. However there are a few modifications which we now mention. 

First note that because a subspace of the stabilizer Hamiltonian contains a degeneracy corresponding to the logical 
qubit and in addition gauge qubits, there can be weight 2 errors acting on the ground state which act non-trivially on 
the degeneracy of the ground state. Thus, whereas before this gave rise to a x-^bo^o term, this can now give rise to an 
0{x^) term which acts trivially on the information in the logical qubit, but which acts nontrivially on the gauge qubit 
degrees of freedom. This means that instead of giving rise to just a projector onto the ground state subspace of order 
x^, the gauge qubits may be split with a strength of 0{x^). Note however, that such interactions, by definition of the 
gauge qubits, act trivially on the information in the logical qubits. For this reason these errors do not give rise to 
new problematic interactions for subsystem codes. Similarly there are now errors which take an excited energy level 
to the ground state, but which do not act on the information encoded into the logical qubits. It is only the errors on 
the logical qubits which matter in deriving the effect of the single qubit perturbations on the encoded information. 



1. Example for the Nine Qubit Subsystem Code 

For the nine qubit subsystem [[9, 1, 4, 3]] code whose stabilizer Hamiltonian is given in Eq. we can carry out an 
analysis of the effect of single qubit perturbations similar to that which we have performed for the five qubit stabilizer 
Hamiltonian. We begin by bounding 69. Here, there are two qubit errors which act non-trivially on the guage qubits. 
However these interactions only act trivially on the information encoded into the logical subsystem of this code. Thus 
we can bound bo in a way similar to the five qubit code: the only interactions which produce an energy shift come 
from the same single qubit interaction acting twice on the subsystem. This implies that |&o| < ^■ 

We can bound 6oi(a) in a matter similar to how we proceeded with the five qubit subspace code. Again the errors 
which are of greatest weight are those which cause the ground state to be excited out of the ground state to the first 
excited subspace. By direct enumeration there are 44 such two qubit errors which contribute to this term. Thus 



max max 6oi ( , a' ) 

i^O a' 



44 

< (34) 



As before we can put this together to bound the effective interaction strength of the the single qubit interactions on 
the encoded quantum information 

11^,11 <|,^,«,^,,(5^)\ (35, 

This becomes less that the original single qubit strength x, when 

7l7 + 60V<l, (36) 
where 7 = 7- This is true for 7 < -^^+^^^^9041 _ 0.0019936. 



C. Single Qubit Errors on Concatenated Distance Three Stabilizer Code Hamiltonians 

We now proceed to analyze the effects of errors on a concatenated stabilizer code Hamiltonian. We will begin by 
considering again the simplest case where the code is a [[n, 1, 3]] code. 

We have shown above that the effect of a single qubit interaction on a stabilizer code of distance d = 3 on the 
ground state is an operator that acts as a logical operator tensored with an operator that acts to excite the ground 
state to a higher energy level. Consider now a concatenated code that has been concatenated only one time. For 
each of the n lowest level stabilizer Hamiltonians, we can apply Eq. (PO)) . If we now focus on the ground state of 
each of these n Hamiltonians (which is 2"-fold degenerate) we see that the effect of our single qubit perturbations are 
now turned into single qubit perturbations on the n encoded qubits tensored with operators that act on the stabilizer 
space of each relevant encoded qubit along with the new stabilizer Hamiltonians that have perturbed ground states. 
But, as we noted when we discussed our system coupled to a bath for a single stabilizer code Hamiltonian, we can 
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still apply our effective Hamiltonian procedure, but now to the encoded qubits. In other words, if the effect of the 
stabilizer code Hamiltonian is to take single qubit interactions of strength x to strength Cx^ on encoded qubits, then 
the effect of these interactions on a concatenated code Hamiltonian for one level of concatenation will act as encoded 
single qubit interactions of strength C(Ca;^)^. 

We now see that our argument allows us to show that single qubit interactions of strength x, at the rth level of 
concatenation produce an error on the encoded qubits of strength less than 

xr < ^{Cxf = X. , (37) 

where a;* = ^ is the threshold coupling strength. Note that the final interaction we obtain will be a logical operator 
tensored with a complicated interaction on the stabilizer spaces of the concatenated code, and, if we started with single 
qubit couplings to a bath, an interaction on the bath. However our argument has always bounded these operators to 
have operator norm 1 and thus Xr represents the true strength of these interactions. Finally note that we also retain 
terms that act on the excited states of the stabilizer Hamiltonians at all levels. However these terms do not couple 
to the ground state. Since we will be interested in the case where the interaction strength to excite to these energy 
levels is exponentially small, these terms will never be relevant if we start our system in the encoded grounded state. 

Thus we see that if the single qubit coupling strength is below x* , then we obtain a doubly exponential suppression 
of these terms in the number of levels of concatenation. Since r levels of concatenation requires n*" qubits, this implies 
that the effective coupling strength shrinks exponentially as function of the number of qubits. 

Finally we note that the new ground state of the concatenated code stabilizer will be close to the ground state of 
the original ground state. At the lowest level of concatenation, our code states will be off by an amplitude of size 
0((-j)^). But these are correctable errors for the next level of concatenation. By the time we get to the rth level of 
concatenation we will have an amplitude which is off by 0(( j)^ ). This will not be correctable, but is exponentially 
small in the number of qubits. There is a corresponding threshold associated with this calculation depending on the 
constant with the uncorrectable component of the amplitude. 

We now point out a few subtleties and drawbacks of our approach. The first important point is to realize that it is 
because we are working a regime where the perturbation series converges that we can move from from level to level in 
our concatenation scheme. If we were to begin to analyze the effect of the perturbation on the entire Hamiltonian on 
n qubits for a quantum concatenated code Hamiltonian, it would not be the case that the perturbation theory would 
converge below some fixed value of the perturbation strength. This is because the requirement that a;||V^|| is smaller 
than half the eignevalue gap would require that x <^ 0{l/n). In other words the perturbation would be guaranteed 
to converge only if the perturbing strength shrunk as a function of the system size. Note that this does not imply 
that such an analysis could not be made to work: it is only a statement that the standard approach to perturbation 
theory would lead to a lack of guaranteed conversion. We get around this problem by analyzing the effect of the 
perturbation on each level of the concatenation and showing how an interaction at one level can be reduced at the 
next level. 

A further more important point is that our analysis shows that the perturbation shrinks as 0(2'') where r is the 
number of levels of concatenation, independent of the code being used. This seems to defy our intuition that instead of 
0(2'") the effect of the perturbation on a [[n, k, d]] code should shrink as O(d^). Note, however, that our analysis, is an 
upper bound on the rate the perturbation shrinks. We obtain 0(2'") in large part because we need to analyze the effect 
of the perturbation with respect to the logical qubits used at the next level of concatenation. In other words, while 
for a [[n, k, d]] the energy splittings to the ground state occur at order d, the fact that there are new eigenstates for 
the perturbed Hamiltonian means that with respect to the unperturbed eigenstates there are nontrivial interactions. 
Whether this is a drawback of our approach or a fundamental obstacle is a question for further investigation. 

In review we have showed that single qubit perturbations on a quantum concatenated code Hamiltonian are expo- 
nentially suppressed as a function of the number of qubits used in the Hamiltonian for single qubit perturbations which 
are weaker than some threshold value. While we have focused on single qubit perturbations, results for two qubit 
perturbations could also be achieved by using a quantum error correcting code designed to correct two qubit errors. 
For the five qubit code, our analysis in Section [ill A 21 gives a threshold value of the coupling of « 0.0020901J and 
for the nine qubit subsystem code, x, « 0.0019936J. 

IV. IMPLEMENTING A CONCATENATED HAMILTONIAN 

Quantum concatenated code Hamiltonians on systems of size n require 0(n)-local interactions. Since such inter- 
actions are not physically reasonable, we must engineer an effective Hamiltonian that yields these interaction terms 
using terms that are physically reasonable, like, for instance two qubit interactions. There are considerable constraints 
on this procedure imposed, for example, by the results of [iOj EU, H^l ■ However the bounds in [13, SU, S^l are not 



13 



robust under procedures which use ancillas to engineer many-qubit interactions using one and two qubit interactions, 
nor are they robust to interactions which are time dependent. Here we focus on the later weakness and show how 
methods with strong pulsed interactions can produce effective evolutions of the desired form [i^, H, EE, El, E3| ■ This 
method is very simplistic and has many draw backs. We include it mostly to show that there is at least some physical 
possibility of implementing quantum concatenated code Hamiltonians. 

A. Bang-bang techniques for many-qubit simulation of stabilizer Hamiltonians 

The basic idea is to use a common technique for engineering many-qubit interaction using strong pulses [1^, IH, E^ ■ 
Consider the quantum circuit given in FIG[TJ This circuit uses two qubit gates conjugated about a single qubit time 
evolution to produce an effective interaction of a four qubit interaction due to the Hamiltonian Heff = X^X^X^X. 
One thing to note about this circuit is that, while we have drawn this as a depth seven circuit, the controlled- nots 
before and after the single qubit gate commute among each other and therefore can be implemented simultaneously. 
In other words physically we can implement this circuit using a circuit of depth of three. A generalization of this 
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FIG. 1: A circuit which is equivalent to the unitary U — exp{—itX ® X ® A (g) A). 

procedure (48l . l49j is easily seen to allow one to implement any Hamiltonian which has a single term corresponding 
to an element of the Pauli group. Depending on how one counts implementation in different single-qubit basis terms, 
a evolution of the form exp(— iPt) where P is a n qubit interaction from the Pauli group can be done with a circuit 
of depth three or depth five. Note that this constant depth is independent of the weight of the Pauli group member 
being implemented. 

Having shown how to obtain an effective interaction of the form exp(— iPt) for P a Pauli group element we 
now turn to how to construct a concatenated stabilizer code Hamiltonian. Note that if we concatenate r times, 
we have rf qubits and the concatenated stabilizer Hamiltonian is the sum of n'' — 1. Stabilizer group elements 
which all commute with each other. Recall that for commuting Hamiltonians, [Hi,H2\ = 0, successive evolutions 
sum: exp{—iHit) ex-p{—iH2t) — exp(— + H2)t). Further note that the concatenated nature of the concatenated 
stabilizer Hamiltonian can be implemented in a highly parallel nature. All of the terms at a given level of the 
concatenation can be implemented in parallel. This implies that using a circuit of depth C(n — l)r where C is 3 or 5 
can be used to implement exp{~iHconcatt)- 

By repeated application of [/ — eyip {—iHconcat^t) we can now produce a system which evolves according to the 
evolution operator U{t) = e:xjp(~iHconcatt) in a stepwise fashion in times of step At. Supposing that our system-bath 
interaction which we are attempting to suppress has a time scale associated with it ottss, then we need to implement 
a single step in this process on a timescale much faster than tss in order for the environment to see an effective 
evolution on the time scale oi tsB- We note that because the depth of our circuit is logarithmic in the number of 
qubits being used, we can expect that for a given strength of system interactions implementing gates we can obtain 
a considerable 

The above method for simulating evolution due to a concatenated stabilizer code Hamiltonian is surely naive. 
It leaves out detailed consideration of how errors in implementing the simulation are dealt with, and a detailed 
understanding of whether the assumption of being able to implement bang-bang operations is valid. This is a subject 
for future research on quantum concatenated code Hamiltonians. 

V. COMPUTING ON QUANTUM CONCATENATED CODE HAMILTONIANS AND 

SELF-CORRECTION 

In this final section we briefly discuss some issues with regards to performing quantum computation on quantum 
concatenated code Hamiltonians. 

One of the main reasons for studying quantum concatenated code Hamiltonians is, that in addition to providing 
a ground state energy manifold which can protect against few qubit perturbations on the ground state, the energy 
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landscape of these systems contains more than just a constant energy gap to the first excited energy state. In Kitaev's 
toric codes in two dimensions [l3| . for example, the ground state is robust to perturbation and there is a constant 
energy gap to the first excited state, but it is possible, by a series of local operations to cause information on the 
ground state to be erred using only a constant amount of energy. In other words, one can excite the system out of 
its ground state at a constant energy, then, at a cost of no energy move across energy levels of the same energy and 
finally relax back down to the ground state while causing an error on the information encoded in the ground state. 
While this does not have a strong effect on the use of such systems as a quantum memory, it has been argued that 
for performing computation on the information encoded in the ground state, such an energy landscape will result in a 
disordering of quantum computation while the computation is being enacted (see [13 )■ In contrast, the four spatial 
dimension toric code has an energy landscape for which an amount of energy proportional to the energy gap is not 
sufficient to cause information encoded into the ground state to be destroyed [ij] . The energy landscape in this later 
case is self- correcting p^[l7l.[25j. 

Here we note that quantum concatenated code Hamiltonians have an energy landscape which is self-correcting. For 
example, for the five qubit concatenated code Hamiltonian discussed in prior section, one can see that the lowest 
weight error that costs the energy gap J in energy to excite out of the ground state is of weight 2. Such an error acts 
on a single lowest block of the code (note that a single qubit error on the ground state for the Hamiltonians we have 
constructed costs rJ energy since it violates all levels of the concatenation.) Having caused an error of weight 2 at 
a cost of energy J, however, an additional two qubit error which brings the system closer to causing an error which 
acts non-trivially on the code costs and energy at least J more (again single qubit errors cost even more.) In other 
words, if we imagine applying small weight errors to the system prepared in its ground state, then an energy barrier 
of height rJ must be cross before the quantum information encoded into the ground state can be acted non-trivially 
upon (r is the number of levels of concatenation in the code.) 

Let us make this statement more rigorous. Suppose that we have a quantum concatenated stabilizer code Hamilto- 
nian using a [[n, 1, d\] code, and this code has been concatenated r times. Consider the set of correctable errors on this 
concatenated code which are corrected using the block nature of the code (that is the lowest level of concatenation is 
corrected first, and then the next level is corrected, etc.) Consider such a correctable error, E, and denote c as the 
number of times in the concatenated code correction where a non-trivial error correction procedure is applied and 
succeeds in correcting the relevant block (disregard corrections which fail or are trivial.) Then E acting on the ground 
state of the code must go to a state of energy at least wJ. This is because each correctable error which is correctly 
fixed must be able to fiip at least one stabilizer generator for the relevant block being error corrected, and this must 
correspond to a minimal change of J for the concatenated quantum code Hamiltonian. Thus we know that if a code 
is concatenated r times, there is an energy barrier of height rJ which must be overcome before an uncorrectable error 
is reached. 

It is interesting to compare the energy barrier for the concatenated stabilizer code Hamiltonian's with the case of 
the four-dimensional toric code[l3|- Here the height of the barrier scales as 0(log A^) for N qubits while for the toric 
code in four-dimensions the barrier of the height scales as 0{N^^^) for N qubits. Whether the slower scaling for the 
concatenated code case leads to less robust behavior is an interesting open question. 

Thus we can imagine how to perform at least some of the basic computations on the quantum concatenated code 
Hamiltonians on the ground state. For example, suppose that we wish to enact an encoded Pauli operator on the 
ground state. Such operators for stabilizer codes are members of the Pauli group. Thus one way to apply this operator 
is to locally apply the relevant Pauli operator to the desired qubits. If one does this in a bang-bang fashion [s^], 
where the strength of this application is much stronger than the natural energy of the quantum concatenated code 
Hamiltonian (0(J)), then one can enact the relevant Pauli group operator on the ground state. However, performing 
bang-bang pulses will inevitably lead to errors during the creation of the local Pauli group operator. For example, 
timing errors may result in over or under rotations of the qubits. For a similar procedure with a simple gapped system 
like Kitaev's two dimensional toric code, such errors will quickly destroy the quantum information in the ground state. 
However, for self-correcting systems, these errors will cause the system to be excited out of the ground state, but if 
the errors are few enough in nature, then the system can quickly relax back toward the ground state and restore the 
quantum information. Note that it is likely that one would still require that the system remain at temperature below 
the energy gap to the first excited state, since occupation of the higher energy levels can cause subsequent quantum 
operations to fail. 

VI. CONCLUSION AND FUTURE DIRECTIONS 

We have introduced a new class of physical systems that provide protection for quantum information encoded into 
the energy levels of these systems called quantum concatenated code Hamiltonians. We have rigorously shown that 
these Hamiltonians offer protection against the effects of perturbing interactions, showing a increase in protection 



15 



which scales exponentially in the number of qubits of the system. Further the energy landscape of these Hamiltonians 
differs from topological quantum systems, which exhibit a energy gap to the first excited states, but can be disordered 
with only a constant amount of further energy once such excited states are created. Quantum concatenated code 
Hamiltonians represent a new method for physically protecting quantum information which is closely connected to 
the traditional method of using concatenation to achieve fault-tolerant quantum computing. This close connection 
should allow for many of the tools and techniques of this standard method to be applied to physical methods for 
building a robust quantum system. 

There is much work, however, that needs to be performed before the ideas described in this paper become a viable 
path toward constructing a large scale robust quantum computer. The first, and most pressing issue, is that while 
we have used one and two qubit interactions gates to construct a bang-bang simulation of the Hamiltonian, it would 
be better if we could find a way to more directly implement these Hamiltonians. One method which is promising 
is to leverage existing results about perturbation theory gadgets [5l|, H^, HI] . However, to date, perturbation theory 
gadgets cannot engineer interactions of great enough simultaneous strength and size to realize concatenated code 
Hamiltonians. 

A second open question is how to perform universal quantum computation on quantum concatenated code Hamil- 
tonians. There are two issues here. One is when the gates to be enacted are performed in a pulsed fashion. While we 
have discussed the enaction of Pauli group operations, a larger set of gates is required in order to achieve universal 
quantum computation. A related question is how the energy landscape of a quantum concatenated code Hamil- 
tonian can be used to suppress computational errors when enacted pulsed gates on the quantum information. A 
detailed understanding of the statistical physics of quantum concatenated code Hamiltonian would be helpful for this 
endeavor [H^. A second approach to performing computation using these Hamiltonians is to use them as building 
blocks for the recently discovered methods for universal adiabatic quantum computation [5l|, HI, [s^, 113, Hll . Such 
adiabatic evolutions may be more suitable for a paradigm where an always on Hamiltonian is used to protect quantum 
information. 

Finally we note that while quantum concatenated code Hamiltonians have Hamiltonians of a highly nonlocal 
nature, there is reason to believe that local few qubit Hamiltonians can give rise to these Hamiltonians. In particular 
we note that concatenated code states fit naturally within the ansatz states of entan glem ent renormalization 59*|. 
Thus concatenated code states serve as candidates, in addition to toric code states [2y], for fixed points in this 
renormalization method. 
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APPENDIX A: KATO'S PERTURBATION THEORY 

Kato ^6^ and Bloch {dT] long ago developed a rigorous perturbation theory that allows one to derive effective 
Hamiltonians for perturbing interactions that is particular useful when the spectrum of the unperturbed Hamiltonian 
has widely separated ei gen values (wide in comparison roughly to the strength of the perturbation.) The reader is 
referred to Kato's work |60j for details of the former perturbation theory and to Bloch, and Jordan and Farhi [6l| 
for details of the later. Here, for completeness, we list the main results of Kato's perturbation theory and in particular 
use the bounds derived in an early version of [53| for errors in truncating this perturbation theory. 

Let Hq be a Hamiltonian with an eigenvalue decomposition Hq = J2i EiHi where Ei the energy of the zth eigenspace 
and Hi is a projector onto this zth, possibly degenerate, eigenspace. We assume that the H^ form a complete set of 
projectors onto the Hilbert space we are working on, H^H^ = (Jy-H^ and J^i^i — ^- We will consider a perturbation 
of Hq by xV, where x is the strength of the perturbation, so that the full Hamiltonian is = Hq + xV. We will 
be concerned with the situation where Hq has discrete energy levels that are well separated from each other, but 
possibly degenerate. In this case as a function of x, the eigenvalues of H split away from their values at Hq in a 
well behaved manner. If we focus on a particular eigenvalue Afc of Hq, then it is possible to obtain an expression for 
the Hamiltonian restricted to the eigenvalues arising from the shift/splitting of this kth eigenvalue of Hq. If the 
entire spectrum of Hq consists of well separated eigenvalues (the exact nature of this condition on the spectrum will 
be described below), then we can obtain an expression for the entire effective Hamiltonian. 
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Begin by defining 



-n, if fc = 



(Al) 



Then the effective Hamiltonian commg from the spHtting of the ith eigenvalue(s) of Hq is 

oo 
m—l 

where 



A, 



(m) 



(m-l) 



(A2) 



(A3) 



and Pi{x) is tfie projector onto tfie eigenvalues of Hx arising from the ith eigenvalue(s) of Ho- The notation ^ 



means that the sum should be taken over all sets of integers kj > such that J2iLi — 
we will often use the first few terms the sum in Eq. (|A2[) . We list these here for ease of reference 
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1. In our calculations 
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Further the projector onto the eigenvalues of arising from the ith eigenvalues(s) of Hq has a perturbative expansion 



P,{x)=Il, + Y,B, 
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where 
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Eqs. (|A2p and (|A5[) are the main expression of importance in this appendix. When is less than half the distance 

from Xi to its nearest other eigenvalue of Hq, Aj, this sum is guaranteed to converge and is correctly reproduces 
the spectrum of H^ arising from the splitting/shifting of the ith eigenvalue of Hq. Sometimes this sum is exactly 
evaluatable and thus one can obtain a closed form expression for H^. On the other hand, the series is more widely 
used to obtain a perturbative expansion in x, truncating the sum in Eq. (jA2p at a fixed m. In this case it is important 
to know the magnitude, and often the form, of the remaining terms in this sum. By an analysis perform in an earlier 
draft of Jordan and Farhi [sl] one can show that if the eigenvalue i is separated from all of other eigenvalues by at 
least A, then if we truncate the series at the pth term, then the remaining terms have an operator norm of no greater 
than 



which, for < 4 is bounded by 



J2 ^""Mm) 
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Similarly for the projector onto the i energy eigenspace we can bound the error as 



< 2 



4a;||V^| 



P+i 



(A9) 
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We note that these bounds are not quite tight, but are reasonable for our purposes. 
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